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Measurement of the azimuthal ordering of charged hadrons with the ATLAS detector 
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This paper presents a study of the possible ordering of charged hadrons in the azimuthal angle relative to the 
beam axis in high-energy proton-proton collisions at the Large Hadron Collider (LHC). A spectral analysis of 
correlations between longitudinal and transverse components of the momentum of the charged hadrons, driven 
by the search for phenomena related to the structure of the QCD field, is performed. Data were recorded with 
the ATLAS detector at centre-of-mass energies of ^ = 900 GeV and ^/s = 1 TeV. The correlations measured 
in a kinematic region dominated by low-pj- particles are not well described by conventional models of hadron 
production. The measured spectra show features consistent with the fragmentation of a QCD string represented 
by a helix-like ordered gluon chain. 



PACS numbers: 13.85.Hd,13.87.Fh 



I. INTRODUCTION 



Inclusive charged-particle distributions have been mea- 
sured in proton-proton {pp) collisions at the LHC for differ- 
ent centre-of-mass energies |ll|-|8|]. These measurements pro- 
vide insight into the strong interaction (QCD) at low energy 
scales and show that predictions of current phenomenological 
models cannot fully describe the measured observables in all 
kinematic regions. As shown in Ref. |;2J, the limitation of 
models is particularly apparent when particles with low trans- 
verse momentum {pr < 500 MeV) are studied. Specifically, 
discrepancies are observed in the description of 

• the charged particle density as a function of pseudora- 
pidity, dA^c/i/d^; 

• the shape of the charged-particle multiplicity distribu- 
tion both at low and high multiplicities; 

• the shape of the charged-particle multiplicity distribu- 
tion for both pt < 500 MeV and pr > 3 GeV; and 

• correlations between the average pj of charged parti- 
cles per event, {pr), and the charged-particle multiplic- 
ity «c/,. 

Some of these discrepancies may be reduced through the 
development of parametrisations for the models of non- 
perturbative QCD and their adjustments (tunes) that better 
match model predictions to the latest measurements with par- 
ticles produced at very low pr- Nevertheless, it is also pos- 
sible that a new formulation of ceitain components of these 
phenomenological models is needed. 

Many of the difficulties in accurately describing observ- 
ables dominated by low-pr QCD phenomena stem from the 
fact that there is often a combination of non-peiturbative ef- 
fects, including soft diffraction, low-y^j- parton scattering and 
hadronisation. These effects act simultaneously in a given 
kinematic region and are difficult to separate experimentally. 
The recent ATLAS study of inclusive charged-particle produc- 
tion ||2] shows that the sensitivity of measurements to some of 



these effects depends on the definition of the kinematic region 
in which the measurements are made. Typically, the more in- 
clusively the sample is defined, or the more soft particle activ- 
ity is added to the measurement, the larger the disagreement 
is between the data and the models. 

Hadronisation, the mechanism of hadron formation from 
quark and gluon cascades, must be taken into account in all 
measurements involving hadrons. The flows of energy, mo- 
mentum and flavour of hadrons approximately follow those of 
partons [9]. The mechanism of hadron formation, however, 
can only be described with the aid of phenomenological mod- 
els. The two main models employed in multi-purpose event 
generators are the stiing (Lund) fragmentation model LI 0.1 and 
the cluster model ul\\ . 



' Full author list given at the end of the article,; 



A proposed modification of the Lund string model yj 
evokes the possibility of the formation of a helix-like struc- 
tured gluon field at the end of the parton cascade. Such a stiTic- 
ture corresponds to the optimal packing of soft gluons in the 
phase space under helicity conservation requirements. Most 
notably, the helix string model imposes coiTelations between 
the adjacent break-up points along a string with observable ef- 
fects in the inclusive pr distribution and in the azimuthal or- 
dering of direct hadrons, i.e. hadrons produced directly from 
the string fragmentation. 

In this paper, the correlations in the charged hadron az- 
imuthal angle (j) lUSll are studied for two different schemes 
of hadron ordering using data recorded with the ATLAS de- 
tector. The results are corrected for detector effects and com- 
pared with the predictions of various Monte Carlo (MC) event 
generators. The paper is organised as follows: the observables 
are desciibed in Section II. Section III describes the ATLAS 
detector and Section IV presents the MC samples used in the 
study. The track and event selection are detailed in Section V. 
The analysis and the uncorrected data distributions measured 
at the center-of-mass energy of y/s=l TeV are presented in 
Section VI. The procedure to correct for detector and selec- 
tion effects and an overview of systematic uncertainties can 
be found in Section VII (additional information is provided in 
the Appendix B). The last section of the paper is devoted to a 
comparison of coiTected data with models and a discussion of 
the results. Appendix A presents results with -^=900 GeV. 



II. OBSERVABLE QUANTITY 



power spectram 



The presence of azimuthal ordering, stemming from the un- 
derlying QCD string structure, is studied for two different 
variants of the helix-like ordered gluon field. As suggested 
in III2I1 . the azimuthal ordering of hadrons should be observ- 
able with the help of a power spectrum defined according to 
the expected structure of the helix field. Assuming the break- 
up of the string occurs via tunneling [ICU, with the partons 
emerging at rest, the azimuthal direction of the hadron's in- 
trinsic transverse momentum coincides with the phase of the 
heUx string in the centre of the string piece which forms the 
hadron. Hence, the azimuthal opening angle of two direct 
hadrons measures the phase difference between two corre- 
sponding points along the string, with the transverse plane 
defined with respect to the string axis. The analysis takes ad- 
vantage of the fact that in soft pp interactions, the QCD strings 
tend to be aligned along the beam axis. 

In close analogy with [12], where the authors assume the 
helix winding is proportional to the rapidity difference be- 
tween hadrons, we define the power spectrum 



"ev event "ch ; 



(1) 



where ^ is a parameter and tJj (0j) is the pseudorapidity (az- 
imuthal angle) of the j-th hadron, A^ev is the number of events, 
and «ch is the number of charged hadrons in the event. The in- 
ner sum runs over charged hadrons in the event and the outer 
sum over events in the sample. 

It is important to note that the form of the helix field is 
not well constrained and that it is possible to find several 
parametrisations of the helix field conforming to the assump- 
tions made in jlZl . One possible scenario 0141] corresponds 
to a static, regular helix structure with the phase difference 
A(p proportional to the amount of internal energy stored in the 
fraction / of string with mass Mq 



A(l) = .y fMo = y kai = ^ ae, 



(2) 



where .5^ is a parameter, K is the string energy density, A^ is 
the difference in the helix phase between two points along the 
string and Al and AE are the length and the energy, respec- 
tively, of the corresponding string piece in its rest frame. The 
energy-distance AE along the string between direct hadrons is 
not directly observable but according to MC studies, the sig- 
nature of the helix-shaped field should be visible with help of 
a very loose approximation of the string by a chain of hadrons 
ordered in pseudorapidity. For the purpose of measuring the 
azimuthal ordering, we thus retain two parameters for each fi- 
nal hadron: the azimuthal angle (j) and the position X along 
the chain, evaluated as 



k<j 



Xj = 0.5 Ej 



L^k, 



(3) 



k=0 



where ^k is the energy of the k-th hadron in the chain, and the 
position of the hadron is associated with the centre of the cor- 
responding string piece. Accordingly, we define an alternative 
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SEi(o)^—j: — ii:exp(/(fi)Xj-0j))|2 
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(4) 



where o is a parameter The inner sum runs over 
pseudorapidity-ordered charged hadrons in the event. 

The presence of a helix-like angular ordering of hadrons of 
either type would manifest itself as a peak in the correspond- 
ing power spectrum; the position of the peak would indicate 
the density of the helix winding. It should be stressed that, 
though formally very similar, Se and Sri are only loosely cor- 
related. A modified form of the helix string implies a differ- 
ence in the experimental signature, such that the presence of 
a heUx gluon field creating a peak in Se does not necessarily 
result in a peak structure in S-q and vice-versa. 

The power spectra can also be expressed as a sum of con- 
tributions from pairs of hadrons 

5,(^) = 1 +-^ £ — £cos(^Arjij-A0,j), 

'^ev event "ch j^j 

SEico) = 1 +-L £ — £cos(a)AXy-A(^ij), (5) 

^ev event "ch j-^j 

where A^ij = 0i — ^j is the opening azimuthal angle between 
hadrons, Arjij = '7i — TJj is their pseudorapidity difference and 
ziXjj =Xi—Xj their energy-distance as defined above. The ab- 
sence of correlations corresponds to 5j)(<^)=l and5£(a))=l. 



III. ATLAS DETECTOR 



The ATLAS detector ||15|] covers almost the entire solid an- 
gle around the collision point with layers of tracking detec- 
tors, calorimeters and muon chambers. It has been designed 
to study a wide range of physics topics at LHC energies. For 
the measurements presented in this paper, the trigger system 
and the tracking devices are of particular importance. 

The ATLAS inner detector has full coverage in (j) and cov- 
ers the pseudorapidity range JT] | < 2.5. It consists of a siHcon 
pixel detector, a silicon strip detector (SCT) and a transition 
radiation tracker (TRT). These detectors are immersed in a 2 T 
axial magnetic field. The pixel, SCT and TRT detectors have 
typical r-(j) position resolutions of 10, 17 and 130 fim, re- 
spectively, and the pixel and SCT detectors have r-z position 
resolutions of 1 15 and 580 /xm, respectively. A track travers- 
ing the full radial extent would typically have 3 silicon pixel 
hits, 8 and more silicon strip hits and more than 30 TRT hits. 

The ATLAS detector has a three-level trigger sys- 
tem: level 1 (LI), level 2 (L2) and the event filter (EF). For 
this measurement, the LI trigger relies on the beam pickup 
timing devices (BPTX) and the minimum bias trigger scintil- 
lators (MBTS). The BPTX are composed of electrostatic beam 
pickups attached to the beam pipe at a distance z = ±175 m 
from the center of the ATLAS detector The MBTS are 
mounted at each end of the detector in front of the liquid- 
argon endcap-calorimeter cryostats at z = ±3.56 m and are 



segmented into eight sectors in azimuth and two rings in pseu- 
dorapidity (2.09 < |tj| < 2.82 and 2.82 < |77| < 3.84). Data 
were taken for this analysis using the single-arm MBTS trig- 
ger, formed from BPTX and MBTS LI trigger signals. The 
MBTS trigger was configured to require one hit above thresh- 
old from either side of the detector. The MBTS trigger effi- 
ciency was studied with a separate pre-scaled LI BPTX trig- 
ger, filtered to obtain inelastic interactions by inner detector 
requirements at L2 and EF [|2l • 



IV. MONTE CARLO SAMPLES 

The analysis has been developed using a PYTHIA 
6.421 lIH] (MC09 tune |[i3]) sample with lai-ge statistics, fully 
simulated fl8j and reconstructed using the standard ATLAS 
reconstruction chain [19]. It should be stressed that the choice 
of the MC tune does not affect the outcome of the analysis 
because the data are corrected in a model independent way. 

For the comparison of the corrected data with the standard 
hadronisation models, MC samples produced with PHOJET 
1.12.1.3512QJ,HERWIG-n- 2.5.1 iHJ aHC-UE7-2/MU900- 
2 tunes lUlXPYTHIA 8.130 (4C tune (H) and a recent tune 
of PYTHIA 6 (AMBT2b |24]) have been used. The MC gen- 
erators PYTHIA and PHOJET employ the Lund string frag- 
mentation model whereas the HERWIG MC is based on the 
cluster model. To study the sensitivity of the power spectra to 
the modification of the string fragmentation model, the data 
are also compared with an alternative implementation of the 
fragmentation process based on the helix string field described 
by Eq. 2111. 



V. DATA SAMPLES 

The measurements reported in this paper were made using 
pp collision data recorded at -Js = 7 TeV. The data were col- 
lected with stable colliding beams at 7 TeV and correspond to 
an integrated luminosity of ~ 190 /ib"' from the beginning 
of the 2010 LHC run 101 ■ A sample of pp collision events 
recorded at y/s = 900 GeV corresponding to an integrated lu- 
minosity of ^ 7 /ib^' [2] was also studied and the results are 
shown in Appendix A. 



The requirements on reconstructed tracks included in the 
analysis are the following: 

• the track is reconstructed by the track reconstruction al- 
gorithm used in [1], with an implicit cut on the trans- 
verse momentum, pr > 100 MeV, and more than 6 hits 
in the silicon detectors; 

• the track has a transverse impact parameter with respect 
to the primary vertex \dQ^\ < 2 mm; 

• the track has a longitudinal impact parameter with re- 
spect to the primary vertex \zq^\ sin(0) < 2 mm; and 

• the track is reconstructed in the pseudorapidity range 
-2.5 <TJ <2.5. 

The requirement on the minimum number of tracks en- 
sures full trigger efficiency fz]- The contributions from the 
beam and non-collision background (cosmic rays and detec- 
tor noise) have been investigated in |2] and found to be neg- 
ligible. Events with multiple primary vertices (less than 0.3% 
of the sample and subsamples defined below) are rejected in 
order to prevent a bias from multiple pp interactions in the 
colliding proton bunches. 



Y.2. Subsample definitions 

The analysis is carried out in parallel on the sample selected 
as described above (henceforth referred to as the "inclusive 
sample") and on two subsamples. 

The first subsample contains events where the transverse 
momentum of any reconstructed track does not exceed 1 GeV 
(max(p7') < 1 GeV). This subsample is called the "low-pr en- 
hanced sample". The effects of parton showering and lateral 
boost are diminished in this selection. The max{pT) < 1 GeV 
requirement selects events with little acollinear jet activity and 
thus the transverse activity is expected to be primarily sensi- 
tive to hadronisation effects. 

The analysis is also performed on a second subsample de- 
fined by a higher track pj cut-off, pr > 500 MeV. This partic- 
ular selection yields a subsample with a significantly reduced 
contribution from diffractive pp interactions (see Table I). This 
subsample is referred to as the "low-pj- depleted sample". 



V.l. Event and track selection 

Events are selected using the following criteria: 

• the event has at least one trigger hit in MBTS; 

• the event has one and only one reconstructed vertex and 
this vertex must have at least three associated tracks; 

• the event has no tracks with pr > 10 GeV; 

• the event has at least six reconstructed tracks («„ > 5) 
passing the requirements below. 



V.3. Selection criteria at particle level 

The comparison between corrected data and MC models re- 
quires an adjustment of the event selection in order to avoid 
a systematic bias. The analysis relies on two main selection 
criteria: the charged-hadron pj- and the charged-hadron mul- 
tiplicity. 

The effect of the low pr cut-off is easily modelled at the 
particle level (MC truth 1261) and the corresponding system- 
atic uncertainty is covered by the uncertainty assigned to the 
correction procedure, described in Section VII. The cut on the 
maximal pr of a track is more selective when applied at the 



particle level as it removes also those MC events which con- 
tain a non-reconstructed high pj- track. The effect is non- 
negligible in the low-pj enhanced sample, where it is cor- 
rected for and a systematic uncertainty is assigned to reflect 
the additional uncertainty. 

The charged-track multiplicity selection criteria on the 
measured data need to be modified in order to emulate MC 
modelling. The adjustment of the charged-track multiplicity 
is done in the following way: for each reconstructed track, a 
random number RND is repeatedly generated according to a 
flat distribution until RND < £tr, £ir being the estimated track 
reconstruction efficiency lIZTIl . The corrected charged-track 
mutiplicity n^/, corresponds to the number of random numbers 
generated for the entire event. The procedure also contains 
an additional correction for the residual content of secondary 
tracks. 

The selection based on n^/, roughly reproduces the average 
charged-particle multiplicity of the particle-level sample se- 
lected with a nf™ cut-off (see Table I). Figure[T]illustrates the 
effect of hadron-level cuts on the true charged-particle mul- 
tiplicity nf™. The choice of the selection cut for the current 
analysis (rich ~ nfl" > 10 ) is aimed at minimising any bias 
in the power spectra related to the loss of events due to the 
detector-level charged-track multiplicity cut-off ritr >5. 

The selection criteria distinguishing between data samples 
are summarised in Table II which also provides the informa- 
tion about the final number of events retained for the analy- 
sis and the mean corrected charged-track multiplicity of each 
sample. The relative fraction of diffractive events (based on 
the nominal cross-section) and the mean charged-particle mul- 
tiplicity of the Monte Carlo selection are given in Table I. 

The final track selection contains 2.8±0.4% of secondary 
tracks according to studies performed on simulated samples. 
Non-primary tracks predominantly arise from hadronic inter- 



xlQ-' 



L. o 1000 



PYTHIAS MC09,\;s = 7 TeV 

— unbiased 



500 




charged-particle multiplicity n 

ch 

FIG. 1 . The impact of the charged-track multiplicity cut-off on the 
true charged-particle multiplicity distribution. The arrow indicates 
the cut on the true charged-particle multiplicity, removing the region 
affected by the loss of events due to the requirement of at least 6 
reconstructed tracks (white area). The shaded area corresponds to the 
detector-level cut n,r > 5. The final analysis selection cut (nc/, > 10) 
is indicated by closed points. 



actions with detector material, photon conversion to electron- 
positron pairs and decays of long lived particles. The aver- 
age reconstruction efficiency for primary charged particles is 
75% for the pr >100 MeV track selection and 84% for the 
Pt >500 MeV selection. The systematic uncertainty due to 
the performance of the track reconstruction is estimated fol- 
lowing the studies performed in Ref. ^. 



Model(tune) 


PYTHIA6(MC09) 


low PT cut 


> 100 MeV 


>500 MeV 


max(p7-) cut 


< 10 GeV 


< IGeV 


< 10 GeV 




Rci, >10 (corrected detector level) 


NT 


31.93 


17.11* 


22.62 


diffractive/total 


3.9±0.1% 


21.4±0.2% 


<0.1% 




n*7 >10 (particle level) 


<" 


31.31 


15.53 


22.26 


diffractive/total 


4.0±0.1% 


22.7±0.2% 


<0.1% 



TABLE 1. Average charged-particle multiplicity N^f" and relative 
fraction of diffractive events for the fully simulated PYTH1A6 
(MC09) MC sample at ^ = 7 TeV. Results are shown for events se- 
lected with the corrected charged-track multiplicity cut-off n^/? >10 
(detector level) and n^™ >10 (particle level) (* indicates before the 
raax(pj) cut correction). 



pp collisions at -^J = 7 TeV, "cA > 10 


PT > 100 MeV 


PT >500 MeV 


max(pr)<10GeV 


max(p7-)<l GeV 


max(pr)<10GeV 


Nev 


Nch 


Nev 


Nch 


Ncr 


Nch 


8,099,211 


34.71 


1,292,389 


17.96* 


4,341,217 


nil 



TABLE II. Number of selected data events and average cor- 
rected charged-track multiplicity, per sample (* indicates before the 
max(/;i7-) cut correction). 



VI. ANALYSIS METHOD 

For the measurement of Se, the selected tracks are ordered 
by pseudorapidity and a pion mass is assigned to each of them. 
According to MC estimates, the charged particle sample con- 
tains about 86% pions, 9.5% kaons, 4% protons/antiprotons 
and a negligible number of leptons (^ 0.5%). The effect of 
assuming a pion mass would need to be taken into account for 
a precision measurement of the position of the signal but its 
impact on the comparison of data with MC models is negligi- 
ble. For the calculation of Sr\ no mass assumption is required. 
The power spectra are measured as {Sr\ ~ 1) and {Se — 1) 
for convenience. 

The uncorrected power spectra are shown in Fig. |2]as a 
function of the azimuthal opening angle (helix phase differ- 
ence) per unit of energy distance (Se) and per unit of pseudo- 
rapidity (S'lj) for the three samples. All angles are expressed 
in radians throughout the paper 
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FIG. 2. The uncorrected data measurements (top: 5^, bottom: 5jj) 
obtained from the data sample collected at y^s = 7 TeV. The mea- 
surement of the inclusive sample is compared to measurements of 
the low-pj- enhanced and low-pj- depleted subsamples. 



decreasing track pr selection range, a feature we may asso- 
ciate with the relative fraction of high pr jets in the sample. 



VII. CORRECTION PROCEDURE AND SYSTEMATIC 
UNCERTAINTIES 

The data are corrected for non-reconstructed charged parti- 
cles with the help of an unfolding technique based on |281 and 
described in the Appendix B. 

The correction for the secondary track content is obtained 
using a random sampling of secondary tracks according to the 
parametrised secondary track rate obtained from fully simu- 
lated MC. The contribution from tracks labeled as secondary 
(^§gsec-j is subtracted from the measured data distribution. Typ- 
ically, it amounts to ^ 6% of the size of the peak in the power 
spectrum. 

An additional correction is applied in the low-pj enhanced 
sample to compensate for the bias introduced by the selection 
cut on max(/?7-) in case the track with highest pj- was not re- 
constructed. Selecting events with exactly one track above the 
threshold, the power spectra (S''^) are calculated using all the 
other tracks in the event. The correction is obtained by sub- 
tracting SP^ from the data in the proportion coiTesponding to 
the probability for the high-pr track being lost in the recon- 
struction. 

All corrections are model independent. The connection pro- 
cedure has been verified by checking the procedure on fully 
simulated MC samples. 



VII. 1. Systematic uncertainties 



MC studies show that the power spectra are sensitive to var- 
ious kinds of coiTelations between particles. The dominant 
peaks seen in both distributions arise due to the jet structure 
and momentum conservation in the hard parton-parton scat- 
tering. The position of these peaks depends on the visible 
energy (in the case of Se) and the pseudorapidity range (in 
the case of 5,j) used in the analysis. Their height is sensi- 
tive to a number of physics processes, notably the amount and 
structure of multiple parton interactions, cross-talk between 
overlapping hadronic systems (colour reconnection), and par- 
ton shower properties. Note, that the values of Se{co = 0) and 
Sri {£, — 0) are identical by definition. They are closely related 
to the average opening angle between particles in the trans- 
verse plane. The presence of azimuthal coiTelations stemming 
from the properties of the gluon field should be visible as an 
additional positive peak or enhancement in the power spec- 
trum and the modification should be more pronounced in the 
low-pj- region, where the fragmentation and parton interac- 
tions have comparable effects on the transverse momentum of 
hadrons. 

The comparison of the uncorrected data obtained in the in- 
clusive event selection and in the low-pr enhanced/depleted 
subsamples shows the size of the peaks diminishing with the 



The principal sources and parametiisations of systematic 
unceitainty associated with the coiTected data are summaiised 
in Table III. The combined systematic uncertainty has the fol- 
lowing components: 

• residual bias of the folding procedure (Appendix B): 
obtained from the comparison of distributions recon- 
structed at the detector level using samples with full 
detector simulation and those obtained with the folding 
technique; 



Source 


Systematic uncertainty [Se{(0), S-qi^)] 


folding procedure 


max(0.003, 3%(5-l)) 


unfolding 


envelope of the residual bias from unfolding 


tracking efficiency 


scaling parameters ±5% 


secondary tracks 


max(0.005, 0.25 ISS'"^]) 


rich cut-off 


variation of the cut-off Wc/i by ± I 


maxipr) < 1 GeV 


SS''^ = (1 - e,r)/£n S''^ , £,r ± 5% 



TABLE III. The parametrisation of the components of the systematic 
uncertainty, per measured point of the corrected power spectrum. 



uncertainty of the unfolding technique: parametrised to 
cover the residual discrepancies in the scaling of 3 fold- 
ing iterations (Appendix B); 

uncertainty on the tracking efficiency estimate: domi- 
nated by the uncertainty on the inner detector material 
description, which translates into a variation of scahng 
factors by 5%; 

uncertainty due to the residual content of secondary 
tracks: set to 25% of the correction applied, with mini- 
mal value of 0.005 (based on MC studies) ; 

uncertainty due to the difference in the charged-particle 
multiplicity selection at the generator level and at the 
detector level: calculated in a model-independent way 
as a variation of the shape corresponding to the change 
of the averaged selected charged-particle multiplicity 
by one unit; and 

the uncertainty in the correction of the bias due the 
max(/:i7-) cut: corresponds to a 5% variation of the track 
reconstruction efficiency. 



All contributions to the systematic uncertainty are combined 
quadratic ally. The negative correlation between track recon- 
struction efficiency and secondary track content is neglected, 
making the uncertainty estimate more conservative. 



VIII. RESULTS 

The results of this analysis obtained for pp collisions at ^/s 
= 7 TeV are presented in this section. Results from this analy- 
sis repeated for -^ = 900 GeV are shown in Appendix A. The 
corrected data are compared with the predictions of several 
commonly used MC models: PYTHIA6, PHOJET, PYTHIAS 
and HERWIG-H-. 

Figure |3] shows the comparison for the inclusive event se- 
lection {rich > 10, pr > 100 MeV and maximal pr < 10 GeV). 
The principal peak structure observed in the power spectra for 
both Se and S-q is roughly reproduced by PYTHIA and PHO- 
JET models and overestimated by HERWIG-H-. The tail of the 
Se distribution around 0.5 < o < 1 rad/GeV is not reproduced 
by any of the models. 
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FIG. 3. Corrected data from the inclusive sample compared to 
particle-level predictions from various MC models using conven- 
tional hadronisation algorithms. The top and bottom plots are for 
the Se and 5(j power spectra, respectively. The error bars correspond 
to the combined statistical and systematic uncertainties. 



FIG. 4. Corrected data from the low-pj- enhanced sample compared 
to particle-level predictions from various MC models using conven- 
tional hadronisation algorithms. The top and bottom plots are for the 
Se and 5tj power spectra, respectively. The error bars correspond to 
the combined statistical and systematic uncertainties. 



Hadronisation effects should become more evident when 
measurements are made in regions of the phase space domi- 
nated by the production of low-pr particles. Figured shows 
the power spectra measured in the low-pr enhanced sample 
{rich > 10, pt >100 MeV and maximal pr < 1 GeV). A 
significant amount of correlations in observed in the data in 
both Se and S-q distributions compared to the PHOJET and 
PYTHIA based models. HERWIG++ gives a seemingly bet- 
ter description for the Se distribution yet it seems its predic- 
tion is more of an artefact of an enhanced single jet structure 
given the fact the model overestimates the measurements in 
the inclusive event selection (Fig.|3]l. The interpretation of this 
measurement in terms of the azimuthal ordering of hadrons re- 
lated to the properties of the gluon field is discussed in Section 
VIII. 1. 

Figure |5] shows the power spectra Se and S-q for the cor- 
rected data and MC predictions in the low-pr depleted region 
(iich > 10, Pt >500 MeV and maximal pr < 10 GeV). In 
principle, this should be the best understood part of the phase 
space, with a suppressed diffractive component, lower sensi- 
tivity to hadronisation effects and best available model tunes. 
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FIG. 5. Corrected data from the low-pj- depleted sample compared 
to particle-level predictions from various MC models using conven- 
tional hadronisation algorithms. The top and bottom plots are for the 
5^ and Svj power spectra, respectively. The error bars correspond to 
the combined statistical and systematic uncertainties. 



However, we find that all models significantly overestimate 
the size of the principal peak structure in both Se and S-q ■ 

Comparison of Figs 3, 4 and 5 show that the azimuthal cor- 
relations are qualitatively different in each subsample, and 
that the standard MC models fail to reproduce these accu- 
rately. Similar conclusion can be drawn for the measurement 
performed at y/s = 900 GeV (Appendix A). 
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FIG. 6. Corrected S-q distributions compared with the particle level 
predictions of PYTHIA6 MC09 for various settings, using the non- 
diffractive pp scattering (full line) as the baseline. Top: low-pj- de- 
pleted subsample. Bottom: low-pj enhanced subsample. The error 
bars correspond to the combined statistical and systematic uncertain- 
ties. 



In the frame of the conventional QCD modelling, we have 
tried to identify the most likely source of the observed discrep- 
ancies. Figure |6] shows the sensitivity of the Sq distribution to 
various components of the QCD modelling implemented in 
PYTHIA 6, taking as a baseline the non-diffractive pp scatter- 
ing scenario (indicated by the full line). In the low-pj- de- 
pleted sample, the size of correlations varies strongly with 
the amount of multiple parton interactions (MPI), of initial 
state radiation (ISR) and with the amount of parton showering. 
The data prefer modelling with enhanced radiation and/or en- 
hanced MPI rate which can be achieved via careful adjustment 
of the relevant model parameters. 

However, such an adjustment typically creates an even 
larger discrepancy in the low-pj- enhanced region, where the 
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parton shower and ISR have smaller influence. Inversely, the 
removal of the MPI and of the diffractive processes increases 
the size of the peak in the modelling of the low-pj- enhanced 
region, but none of these rather extreme variations lead to 
a good agreement with the measured Srj distribution, while 
creating a huge discrepancy in the low-pj- enhanced region. 
We conclude that for both measured spectra, the conventional 
models fail to describe the low-pj- enhanced region, where the 
data are consistenly showing a larger and broader peak struc- 
ture. 



VIII. 1. Alternative fragmentation model 

The question of compatibility of the measurements with 
the azimuthal ordering signal originating from the underly- 
ing structure of the QCD field is studied using the PYTHIA6- 
based helix string model implementation Il25ll of the modi- 
fied helix string scenario corresponding to the Se definition 
(Eq. 2). The comparison of corrected data with the PHOJET 



modelling of pp interactions, interfaced alternatively with the 
standard string fragmentation and with the helix string frag- 
mentation, is shown in Figs |7] and [H] It is seen that the helix- 
like gluon ordering improves the description of the data in the 
inclusive sample for Se as it generates higher values in the re- 
gion of O G (0.5, 1). In the low-pj- enhanced sample the data 
are more strongly peaked than the model. This indicates the 
data have a more jet-like structure than predicted by PHOJET 
which means the model may need readjustment beyond the 
fragmentation part. 

It is possible that the original helix string proposal I112ll pro- 
vides an improved description of the Srj measurement but we 
cannot verify this hypothesis due to the absence of the MC 
implementation of this variant of the helix string model for 
parton configurations containing hard gluons. 

In principle, the position of the azimuthal ordering signal 
allows a measurement of the characteristic density of the he- 
lix string winding (parameter y in Eq. 2 ). The data indicate 
y in the range 0.5 to 0.9 rad/GeV. However, the interplay be- 
tween string fragmentation and the other processes (such as 
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FIG. 7. Corrected 5^ (top) and 5tj (bottom) distributions com- 
pared with the particle level predictions of PHOJET, interfaced to the 
PYTHIA6-based alternative fragmentation using helix string model 
I25l1 (setting .5^= 0.7 rad/GeV, see Eq. 2), in the inclusive event se- 
lection. The erTor bars correspond to the combined statistical and 
systematic uncertainties. 



FIG. 8. Cort'ected 5^ (top) and S-q (bottom) distributions com- 
pared with the particle level predictions of PHOJET, interfaced to the 
PYTHIA6-based alternative fragmentation using helix string model 
125 1 (setting S^= 0.7 rad/GeV, see Eq. 2), in the low-p/- enhanced 
subsample. The error bars correspond to the combined statistical and 
systematic uncertainties. 
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parton showering, colour reconnection and multiple interac- 
tions) is not sufficiently understood to allow a more precise 
estimate (the measurement has to rely on MC description for 
the correction of the longitudinal boost). 



IX. CONCLUSIONS 

A measurement of the ordering of charged hadrons in the 
azimuthal angle with the ATLAS data recorded from proton- 
proton collisions at -/s = 7 TeV and ^/s = 900 GeV has been 
presented. 

A spectral analysis of correlations between the opening az- 
imuthal angle and the longitudinal separation of the charged 
hadrons was performed by measuring the Se and Sri power 
spectra. These measurements were done in three kinematic re- 
gions (inclusive, low-pj enhanced and low-pj- depleted sam- 
ples) that were specifically defined to help assess the potential 
contribution of hadronisation effects to the power spectra by 
varying the levels of competition between hadronisation and 
other QCD effects. 

The results were compared with the expectations of various 
MC event generators. Predictions generated by the MC mod- 
els employing the standard Lund string fragmentation model 
roughly reproduce the data in the inclusive sample. 

The models systematically overestimate the size of corre- 
lations in the low-pr depleted sample, where the observables 
are sensitive to the multiple jet structure of events (due to the 
presence of underlying event and/or parton showering). 

For observables measured in the low-pr enhanced sam- 
ple, none of the models investigated describes the data ade- 
quately. A study showing the impact of extreme variations in 
MC model parameters that are known to contribute to soft- 
QCD effects demonstrates that, although some improvement 
of predictions for Se and 5,j can be achieved in the low-pr 
enhanced sample, it is still far from satisfactory. 

The measurement of Se in the kinematic region dominated 
by low-pr particles shows features similar to those seen in 
models in which the fragmenting QCD strings are represented 
by helix-like ordered gluon chains. These measurements sug- 
gest that the inclusion of such azimuthally-ordered fragmen- 
tation effects could be one factor in improving current models 
of soft particle production and hadronisation. 
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APPENDIX A 

The analysis is repeated for ATLAS data collected at ^/s = 
900 GeV. The number of selected events and the average 
corrected charged-track multiplicity are shown in Table IV. 
For comparison, the average charged-particle multiplicity of 
PYTHIA6 MC09 samples is shown in Table V, together with 
the fraction of diffractive events in the sample (the estimate 
is based on the nominal cross section obtained from the gen- 
erator). The corrected data are compared to the particle level 
prediction of various models in Fig. |9] 



pp collisions at ^ = 900 GeV, rich >10 


PT >100MeV 


PT >500 MeV 


max(/77-)<10GeV 


max(p7-)<l GeV 


max(p7-)<10GeV 


Nev 


K-h 


A'.. 


Nch 


Nev 


Nch 


224717 


24.8 


59880 


17.7* 


68456 


16.3 



TABLE IV. Number of selected events and average corrected 
charged-track multiplicity, per sample (* indicates before the 
max(p7-) cut correction). 
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APPENDIX B 

The deconvolution technique employed in this analysis is 
a model-independent procedure suitable for observables with 
linear dependence on the track reconstruction efficiency. The 
dependence of the shape of the measured distribution on the 
fraction of reconstructed tracks is studied by convolution with 
the track reconstruction efficiency matrix ("double-folding"): 
the reconstructed tracks are randomly rejected from the sam- 
ple according to the parametrised reconstruction efficiency 
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FIG. 9. Corrected data compared to particle level predictions from various MC models using conventional hadronisation algorithms. The left 
and right plots are for the 5^ and 5,j power spectra, respectively. Top: inclusive event selection (pj >100 MeV, max(p7-) <10 GeV). Middle: 
event selection with enhanced soft component (pj >100 MeV, max(p7-) <1 GeV). Bottom: event selection with reduced soft component 
{pT >500 MeV, maxipr) <10 GeV). ATLAS data collected at y^ = 900 GeV and selected with n^i, > 10. The error bars correspond to the 
combined statistical and systematic uncertainties. 
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Model(tune) 


PYTH1A6(MC09) 


low pj cut 


>100MeV 


>500 MeV 


max(pT) cut 


< 10 GeV 


< IGeV 


< 10 GeV 




n^/, >10 (corrected detector level) 


NfT 


21.91 


15.28* 


15.2 


diffractive/total 


3.8(±0.1)% 


14.7(±0.2)% 


<0.01% 




"fr >10 (particle level) 


<" 


21.45 


15.47 


15.75 


diffractive/total 


3.7(±0.1)% 


13.8(±0.2)% 


<0.01% 



TABLE V. Average charged-particle multiplicity 



ch 



and rela- 



tive fraction of diffractive events for the fully simulated PYTH1A6 
(MC09) MC sample at 0^ = 900 GeV. Results are shown for 
events selected with the corrected charged-track multiplicity cut-off 
rich >10 (detector level) and n^^i^ >10 (particle level) (* indicates 
before the max(p7-) cut correction). 



and the power spectra are recalculated. The procedure is re- 
peated two more times, so that a sequence of three folding iter- 
ations is available for each measured distribution. The folding 
iterations, together with the measured distribution and the de- 
convoluted (MC truth) distribution, obey simple scaling rules: 

• the size of correlations (along the S-1 axis) scales lin- 
early with the multiplicity (physically it depends on the 
fraction of correlated pairs which have quadratic mul- 
tiplicity dependence, partially compensated by the nor- 
malisation factor in Eqs. 1 and 4). 

• the shape of the power spectrum Se scales Unearly in 
(0 with the fraction of energy removed from the hadron 
chain (in S-q, the position of the peak stays nearly con- 
stant). 

A typical example of the variation of the power spectrum with 
the number of applied folding iterations is shown in Fig.llOlon 
a PYTHIA6 sample with inclusive event selection. For better 



illustration of the scaling symmetry, secondary tracks are re- 
moved from the sample using the generator level information. 
The deconvolution of the power spectrum requires finding a 
set of scaling factors /,. , /^ (x = O , (^ ) which fuUfil the require- 
ment 

5,_i(x)-l = -^[5,(y)-l], 0<x<10, 0</<4. (Bl) 

where St designates the i-th folding iteration of the power 
spectrum (i=0 for uncorrected data, -1 for corrected data). 

The origin of observed correlations may be diverse, with 
variable dependence on the track reconstruction efficiency. A 
better precision can therefore be obtained by splitting the ob- 
served distribution empirically into components and by esti- 
mating the scaling factors per component. 

The power spectra calculated in this paper can, to good ap- 
proximation, be fit with a combination of a Landau distribu- 
tion 1I29I1 (to describe the resonant peak structure), an expo- 
nentially falling background (driven by local momentum con- 
servation of adjacent hadron pairs) and a polynomial of first 
or second degree to describe the upper tail of the distribution. 
In practice, the subtraction of a single background term and 
determination of 2-3 effective scaling factors is sufficient to 
describe the difference between a pair of folding iterations; 
the replacement of the fx with an attenuated scaling factor 
fx^fx = {fx+x)/{\+x) helps smooth the transition into 
the upper part of the spectrum. 

The knowledge of a single pair of consecutive iterations is 
sufficient to calculate the scaling factors. Typically, they are 
determined using the difference between uncorrected data and 
the 1st folding iteration (1-^0). The stability of the scaling 
factors is verified by applying them on the folding iterations 
of higher orders (2— >1,3— ?>2) while the residual discrepancies 
serve as a basis for the estimate of the associated systematic 
uncertainty (Fig. fTTT i. The deconvoluted power spectrum is 
obtained with the help of Eq. Bl as the extrapolated (0— >-l) 
member of the sequence of folding iterations. 
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FIG. 10. Power spectrum calculated at the generator level, detector 
level and following three folding iterations as described in the text. 
Secondary tracks are removed from the samples in Figs 10-11. 



FIG. 1 1 . The closure test of the scaling deconvolution method. 
Follow-up of the scaling precision over several folding iterations al- 
lows to verify the stability and to estimate the associated systematic 
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[28] J. W. Monk and C. Oropeza-Barrera, arXiv: 1111 .4896v2 

[NIM (to be published)]. 
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